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Within the framework of an SO(IO) GUT model that can accommodate both the LMA and LOW 
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Over the last few years the evidence for neutrino oscil¬ 
lations between the three known neutrino flavors {ue, I'/r, 
and Vr) has become increasingly convincing. The at¬ 
mospheric neutrino flux measurements from the Super- 
Kamiokande (Super-K) experiment exhibit a deficit of 
muon neutrinos which varies with zenith angle (and 
hence baseline) in a way consistent with —> Vr oscil¬ 
lations [Q. In addition, recent combined evidence from 
Super-K and the SNO experiments 1^ indicate that some 
electron-neutrinos from the sun are oscillating into muon 
and/or tau neutrinos. While the atmospheric neutrino 
data with oscillations points to a small region 

of the mixing parameter space [Q, the solar neutrino 
data is consistent with at least two regions of parame¬ 
ter space 1^, corresponding to either the Large Mixing 
Angle (LMA) or to the LOW MSW 0 solution. 

Neutrino oscillation data constrain Grand Unified The¬ 
ories (GUTs) which provide a theory of flavor and relate 
lepton masses and mixings to quark masses and mix¬ 
ings. It is known that the presently implied neutrino 
mass scales can be accommodated naturally within the 
framework of GUTs by the seesaw mechanism [Q. In 
practice finding an explicit GUT model for the LMA so¬ 
lution has been found challenging. However, one exam¬ 
ple has been constructed by Barr and one us |^. In this 
model the Dirac and Majorana neutrino mass matrices 
are intimately related. It has been shown by the present 
authors that by varying the Majorana mass matrix 
parameters any point in the presently-allowed LMA re¬ 
gion can be accommodated. In the present paper we show 
that the LOW region can also be realized in the model by 
choosing an appropriate texture for the right-handed Ma¬ 
jorana mass matrix. In addition to the Majorana mass 
matrix, we also spell out the Dirac mass matrices, list 
the values of the associated input parameters and give 
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results for the quark and charged lepton sectors. For each 
choice of the Majorana matrix we discuss the oscillation 
predictions. We use our results to illustrate how Neu¬ 
trino Superbeams and Neutrino Factories can further 
test this GUT model. Our results suggest that, indepen¬ 
dent of which is the preferred solution (LMA or LOW), 
Neutrino Superbeams and Factories will be necessary to 
identify the correct model, and that further investment 
in developing them should be encouraged. 

Within the framework of three-flavor mixing, the fla¬ 
vor eigenstates Va (a = e, r) are related to the mass 
eigenstates Vj {j = 1,2,3) in vacuum by 

Va = X ^ = Umns^m ( 1 ) 

j 

where U is the unitary 3x3 Maki-Nakagawa-Sakata 
(MNS) mixing matrix Q times a diagonal phase matrix 
= diag(e®^U 1). The MNS matrix is conven¬ 

tionally parametrized by 3 mixing angles (6*23,6*12,013) 
and a GP-violating phase, 5cp' 

Umns = 

C12C13 S12C13 S13U \ 

— S12C23 — Ci 2S23S13^ C12C23 — Si2S23Sl3^ S23C13 j 

S12S23 — Ci2C23Sl3^ —C12S23 — Si2C23Sl3^ C23C13 / 

. (2) 

where cjk = cosOjk, Sjk = sinOjk and ^ 

The three angles can be restricted to the first quadrant, 
0 < dij < 7112, with 5cp in the range —tt < 5cp < tt, 
though it proves advantageous to consider @13 in the 
fourth quadrant for the LMA solutions. 

The atmospheric neutrino oscillation data indicate 
that 

Aruon — 3.0 X 10“^ eV^, 

(3) 

sm^26atm = 1-0, (> 0.89 at 90% C.I.), 

where Amf^ = mf — rrVj and mi, m2 and m3 are 
the mass eigenstates. The atmospheric neutrino oscilla¬ 
tion amplitude can be expressed solely in terms of the 


B,] bri ght@fn a.]. gov 
ggeer@fnal.gov| 









2 


Umns matrix elements and is given by sin^ ‘ 29 atm = 
4|C/^3p(l - \U^3\^) ~ 4|17^3p|C/r3p = 4sin2 2023. The 
approximation is valid because |17e3| is known to be small 


The solar neutrino oscillation data from Super-K in¬ 
dicate that, for the LMA solution, the allowed region is 
approximately bounded by 

- (2.2 - 17) X eV^ 
sm^29,oi ^ (0.6-0.9), 
while for the LOW solution, 

Amii ~ (0.3-2) X 10-'^ eV^, 

tan2 6»i2 (0.6-1.2), 


(4) 


(5) 


where the solar neutrino oscillation amplitude is 
sin^ 29soi = 4|[/ei|2(l - |C/,i|2) ~ 4ll7eip|17e2p, while 

tan2 6»i2 = \Ue2/Uel\‘^. 

The GUT model we consider is based on an S'O(IO) 
GUT with a [/(I) X Z 2 X Z 2 flavor symmetry. The model 
involves a minimum set of Higgs fields which solves the 
doublet-triplet splitting problem. The Higgs superpoten¬ 
tial exhibits the U(1) x Z 2 x Z 2 symmetry which is used 
for the flavor symmetry of the GUT model. Details of 
the model can be found in We simply note that the 
Dirac mass matrices U, D, N, L for the up quarks, down 
quarks, neutrinos and charged leptons, respectively, are 
found for tan /3 ~ 5 to be 


/V 

0 

0 \ 


d 

5’e^* \ 

U = 0 

0 

e/3 

, D= 6 

0 

a + e/3 1 

Vo 

-e/3 1 / 

\5'e^* -e/3 

1 ) 


/V 

0 

0 \ 


6 

0 

0 


, L={ S 

0 

Vo 

e 

1 / 

V 

cr -1- e 



N = 


( 6 ) 

where U and N are scaled by Mu, and D and L are scaled 
hy Md- All nine quark and charged lepton masses, plus 
the three GKM angles and GP phase, are well-fitted with 
the eight input parameters 


Mu ^113 GeV, 
cr = 1.78, 

S = 0.0086, 

(j) = 126°, 


Md ^ 1 GeV, 
e =0.145, 

V = 0.0079, 
r] = 8x 10-6, 


(7) 


defined at the GUT scale to fit the low scale observables 
after evolution downward from Kqut'- 


mt(?7it) = 165 GeV, 
TO„(1 GeV) = 4.5 MeV, 

Vus = 0.220, 

Vcb = 0.0395, 


mr = 1.777 GeV, 
= 105.7 MeV, 
nie = 0.511 MeV, 
Scp = 64°. 


( 8 ) 


These lead to the following predictions: 

mbimb) = 4.25 GeV, mcijnc) = 1.23 GeV, 
m«(l GeV) = 148 MeV, md{l MeV) = 7.9 MeV, 
|K&/Kh| = 0.080, sin2/3 = 0.64. 

(9) 


With no extra phases present, the vertex of the GKM 
unitary triangle occurs near the center of the presently 
allowed region with sin 2d ~ 0.64, comparing favor¬ 
ably with recent results The Hermitian matrices 

D^D, and VtV are diagonalized with small left- 
handed rotations, Uu, Up, Un, respectively, while L^L 
is diagonalized by a large left-handed rotation, Up- This 
accounts for the small value of |Uc6| = |(4?7p)cb|, while 
IG^sl = will turn out to be large for any rea- 

sonable right-handed Majorana mass matrix, Mr |Q. 

The effective light neutrino mass matrix, M,j, is ob¬ 
tained from the seesaw mechanism once M^ is specified. 
While the large atmospheric neutrino mixing <-> v-r 
arises primarily from the structure of the charged lepton 
mass matrix, the structure of Mr determines the type of 
i^T solar neutrino mixing. 

To obtain the LMA solution requires some fine-tuning 
and a hierarchical structure for Mr, but this can be ex¬ 
plained in terms of Froggatt-Nielsen diagrams . Here 
we restrict our attention to a slightly less general form 
for Mr than that considered in Q and 0: 

( b^rj^ —berj arj \ 

-beri -e I Ar, (10) 

arj -e 1 / 

where the parameters e and 77 are those introduced in 
Eq.(^ for the Dirac sector. This structure for Mr can 
be understood as arising from one Higgs singlet which 
induces a AL = 2 transition and contributes to all nine 
matrix elements while, by virtue of its flavor charge as¬ 
signment, a second Higgs singlet breaks lepton number 
but modifies only the 13 and 31 elements of Mr. As 
shown in detail in Q, we can introduce additional GP 
violation by assigning a relative phase to the two lepton 
number breaking Higgs singlets, whereby we set 

a = b — a . (11) 

On the other hand, we find the LOW solution can be 
obtained with the simple hierarchical structure for Mr, 

/ e d 0 \ 

MR=id 0 0 Afl, (12) 

Vo 0 1/ 

where by the flavor charge assignments, one Higgs singlet 
inducing a AL = 2 transition contributes to the 12, 21 
and 33 elements, while a second Higgs singlet also breaks 
lepton number but contributes only to the 11 matrix el¬ 
ement. For simplicity we keep both d and e real, since 
the leptonic GP phase is inaccessible to measurement for 
Am 2 i values in the LOW region. 

For either the LMA or LOW version, M^, is then ob¬ 
tained by the seesaw formula [^, M,, = N'^M^^N. With 
Ml, complex symmetric, both M^M^, and M^, itself can 
be diagonalized by the same unitary transformation, C/^, 
where in the latter case we find 

U'^Mi.Uv = diag(TOi, -m 2 , m 3 ). 


(13) 
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FIG. 1: The viable region of GUT parameter space consistent with the present bounds on the LMA MSW solution. Contours 
of constant sin^ 20 i 2 are shown together with (a) contours of constant and (b) contours of sin^ 2 ^ 13 . 


TABLE I: List of four points selected in the LMA allowed parameter region to illustrate the neutrino oscillation param¬ 
eter predictions of the GUT model. Here the CP phase Scp arises from in L alone, as no phase (p' has been introduced in Mr. 


Am'' 


''21 (eV^) 


Ami 


tarn 9i2 


sin^ 2012 


sin 2023 


sin"^ 2013 


Sci 


1.0 

2.0 

6.5 X 10“® 

3.2 X 10“® 

0.49 

0.88 

0.994 

0.0008 

-4° 

1.2 

2.8 

3.3 X 10“® 

3.2 X 10“® 

0.43 

0.84 

0.980 

0.0038 

- 1 ° 

1.6 

2.9 

6.1 X 10"® 

3.2 X 10“® 

0.35 

0.77 

0.998 

0.0015 

-3° 

1.7 

2.7 

10.9 X 10“® 

3.2 X 10“® 

0.32 

0.73 

0.996 

0.00008 

-14° 

1.7 

3.4 

4.0 X 10“® 

3.2 X 10“® 

0.33 

0.75 

0.992 

0.0033 

- 2 ° 

2.2 

3.5 

8.8 X 10“® 

3.2 X 10“® 

0.24 

0.63 

0.996 

0.0008 

-4° 

real light neutrino masses, Ui, 

can not be arbitrarily 

in the LOW version. 

In either version, the third param- 


phase transformed and is uniquely specified up to sign 
changes on its column eigenvectors 0. Hence Umns 
is found by applying phase transformations on Uj^Ui, to 
bring Uj^Ui, into the parametric form of Eq. whereby 
the el, e2, ^3 and t3 elements are real and positive, the 
real parts of the fj.2 and rl elements are positive, while 
the real parts of the /rl and r2 elements are negative. The 
inverse phase transformation of that applied on the right 
can then be identified with the Majorana phase matrix, 
of Eq. (||). The evolution of the predicted values 
between the GUT scale and the low scales can be safely 
ignored Q, since tan/3 ~ 5 is moderately low and the 
neutrino mass spectrum is hierarchical with the opposite 
CP parities present in Eq. dl)- 

We can now examine the viable region of GUT 
model parameter space that is consistent with either 
the LMA or LOW solar neutrino solution, and ex¬ 
plore the predicted relationships among the observables 
sin^2d23, sin^ 20 i 2 , sin^ 2 di 3 , Scp, Am| 2 , and Am^i. 
We shall emphasize here the simpler cases in which there 
are, in effect, only two additional real dimensionless GUT 
model parameters, a and b in the LMA version or d and e 


eter An sets the scale of Am'^ 2 - The more general CP 
results obtained for the LMA solution with the presence 
of a complex parameter a in Eq. ( 0 ) have been explored 
in detail in 0 . 

The viable region of GUT model parameter space con¬ 
sistent with the LMA solar solution is shown in Fig. H. 
Both parameters a and b are constrained by the data to 
be close to unity, with 1.0 5 a 2.4 and 1.8 P: b < 5.2. 
Superimposed on the allowed region. Fig. 0 (a) shows 
contours of constant sin^ 2 di 2 and contours of constant 
Amf^. Figure 0(b) similarly displays the allowed region 
with contours of constant sin^ 2 di 2 and sin^ 2^13 super¬ 
imposed. The nearly parallel nature of the contours of 
Am^i in (a) and sin^ 2 di 3 in (b) indicates a strong corre¬ 
lation between them. As the predicted Am^i increases, 
the predicted sin^ 2 di 3 decreases. Note that if the LMA 
solution is indeed the correct solution, KamLAND is 
expected to provide measurements of Am^i and sin^ 0 i 2 
to a precision of about 10% From these measure¬ 
ments the model parameters a and b can be determined 
from Fig. 0(a). Figure 0(b) can then be used to give a 
prediction for sin^ 20 i 3 with a precision also of order 10 %. 
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FIG. 2: The viable region of GUT parameter space consistent with the present bounds on the LOW MSW solution for (a) 
negative d and (b) positive d. Contours of constant sin^ 2613 , sin^ 2612 , and Am^i are shown. 

TABLE IL List of four points selected in the LOW allowed parameter region to illustrate the neutrino oscillation parameter 
predictions of the GUT model. 


d 

e 

Amii (eU®) 

Amia (eU®) 

tan® 012 

sin® 2012 

sin® 2023 

sin® 2013 

- 4.2 X 10“® 

10.0 X 10“® 

1.20 X 10"’’ 

3.0 X 10"® 

0.56 

0.906 

0.911 

0.028 

-4.1 X 10“® 

4.0 X 10"® 

0.52 X 10"’’ 

3.0 X 10"® 

0.81 

0.975 

0.899 

0.027 

-3.6 X 10“® 

3.0 X 10"® 

0.64 X 10"’’ 

3.0 X 10"® 

0.86 

0.980 

0.898 

0.030 

3.6 X 10~® 

5.0 X 10"® 

0.98 X 10"’’ 

3.0 X 10"® 

1.00 

0.999 

0.914 

0.0016 

5.3 X 10“® 

10.0 X 10"® 

0.50 X 10"’’ 

3.0 X 10"® 

0.82 

0.989 

0.912 

0.0039 

5.0 X 10~® 

13.0 X 10"® 

0.85 X 10"’’ 

3.0 X 10"® 

0.70 

0.966 

0.918 

0.0033 


In Table | we have selected six points in the LMA al¬ 
lowed parameter region to illustrate the neutrino oscil¬ 
lation predictions of the GUT model. The correlations 
noted above are evident. It is also striking how nearly 
maximal are the values for the atmospheric mixing pa¬ 
rameter, sin^ 2023- This apparently arises not from some 
additionally imposed symmetry but rather from the fine 
tuning between the right-handed Majorana and Dirac 
neutrino mass matrices, cf. Eqs. and How¬ 

ever, if an additional phase is incorporated into Mu for 
this LMA case as indicated in Eq. (|^), the maximal- 
ity of the atmospheric mixing is decreased to the lower 
bound in Eq. as |(5cp| approaches 50°. See Q for 
more details. 

Turning now to the GUT model version for the LOW 
solution, we find that there are two parametric regions 
shown in Fig. 2 for the presently allowed solutions, cor¬ 
responding to —4.8 d X 10® —2.2, 0 $ e x 10® ^ 13 

and 2.0 ^ d x 10® < 6.0, 0 < e x 10® < 25. Here no dra¬ 
matic correlation between sin® 20 i 3 and Am^i exists, so 
we have plotted contours of Am^i, sin® 26*12 and sin® 26*13 
on the same figures. If KamLAND fails to see a signal for 
the LMA region while Borexino 0 , for example, iden¬ 
tifies oscillations corresponding to the LOW region and 
determines sin® 29i2 and Am^i with nearly 10 % precision 
Ip^ , sin® 26*13 will be specified up to a two-fold ambigu¬ 


ity in the GUT model in question. A first measurement 
of sin® 2013 would resolve the ambiguity, and a precise 
measurement would test the model. For the negative d 
version, a SuperBeam facility capable of probing down to 
sin® 2013 — 0.003 will be able to test the model, while for 
the positive d version the complete parameter space can 
only be tested with a Neutrino Factory. Table H gives 
the relevant mixing solutions for a set of six points. In 
contrast to the LMA results with small CP phases, we 
see that the atmospheric mixing for the LOW solution is 
large but not nearly so maximal. 

In conclusion, we have studied predictions for a partic¬ 
ular but representative GUT model that can accommo¬ 
date both the LMA and LOW solar neutrino solutions. 
We find that precise measurements of sin® 20 i 2 , Ato^i, 
and sin® 20 i 3 are needed to test the theory. Given the 
observed near maximal value of sin® 2023 the LMA solu¬ 
tion, which requires some fine tuning of the M}i matrix, 
is favored by the model. The model then predicts that 
the CP phase 6 cp is small and sin® 20 i 3 < 0.006. For the 
LOW solution which requires no fine tuning, sin® 20 i 3 
can be as small as this, or an order of magnitude larger, 
depending upon the sign of the d model parameter in 
Mji. Our work suggests progress on testing GUTs can be 
made with Neutrino Superbeams, but ultimately a Neu¬ 
trino Factory will be needed to help identify the correct 







5 


model. 
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